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Abstract 



A number of recent experiments have showed that surfactants can modify 
the growth mode of an epitaxial film, suppressing islanding and promoting 
layer-by-layer growth. Here a set of coupled equations are introduced to de- 
scribe the coupling between a growing interface and a thin surfactant layer 
deposited on the top of the nonequilibrium surface . The equations are derived 
using the main experimentally backed characteristics of the system and basic 
symmetry principles. The system is studied using a dynamic-renormalization- 
group scheme, which provides scaling relations between the roughness expo- 
nents. It is found that the surfactant may drive the system into a novel phase, 

in which the surface roughness is negative, corresponding to a flat surface. 
PACS numbers: 68.55 -a, 68.35 Fx, 64.60 Ht 
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I. INTRODUCTION 



Lately there is much theoretical interest in the statistical properties of nonequilibrium 
interfaces. Most of the growing interfaces naturally evolve into self-affine structures; the 
surface morphology and the dynamics of roughening exhibit simple scaling behavior despite 
the complicated nature of the growth process ||T|-0 . In particular, much attention has focused 
on different models to describe thin-film growth by molecular-beam epitaxy (MBE) 0-0. 

Under ideal MBE conditions the primary relaxation mechanism is surface diffusion, which 
conserves the mass of the film. Experimentally both lattice strain and surface free energy de- 
termine whether the film undergoes layer-by-layer growth, islanding, or layer-by-layer growth 
followed by islanding. In experiments involving growth of Ge on Si (100) surface layer-by- 
layer growth is limited to 3-4 monolayers (ML) due to the lattice mismatch between Si and 
Ge and is followed by formation of unstrained Ge islands. It was shown recently that is- 
landing in the Ge/Si system can be suppressed effectively by use of a surfactant monolayer, 
changing the growth mode from island growth to layer- by-layer growth . Suitable surfac- 
tants such as As and Sb strongly reduce the surface free energy of both Si and Ge surfaces 
and segregate at the surface during growth. 

In this paper we study the generic problem of nonequilibrium roughening of an interface 
covered by a thin surfactant layer (see Fig p. Building on experimental results and general 
symmetry principles, a set of nonequilibrium equations are proposed to describe the growth 
of an interface coupled to the fluctuations in the surfactant coverage. The analytic study 
of these equations indicates that the surfactant changes drastically the morphology of the 
interface in 2+1 dimensions. In particular, the coupled system supports the existence of a 
novel phase characterized by negative roughness exponent, which can be identified with a 
morphologically fiat surface. A summary of the main results were presented in an earlier 
publication |16l . 



The paper is organized as follows: Section II presents a short review on the experimental 
studies. Section III introduces the key elements of the proposed nonequilibrium theory. 
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Section IV presents the main results of the dynamic renormahzation group (DRG) analysis, 
and finally Section V refiects on open problems and possible future developments. 

II. EXPERIMENTAL STUDIES ON SURFACTANT MEDIATED GROWTH 

The ultimate goal in crystal growth by MBE or other vapor-phase techniques is to 
control and influence the growth mode of a thin fllm. At high temperatures, necessary 
to obtain sufficient mobility of the surface atoms, the growth mode is determined by the 
interface and surface free energies and the lattice strain. Lattice strain is specially relevant 
in heteroepitaxial growth, when attempt is made to combine different type of atoms in a 
layered structure. 

Depending on the surface, interface and heteroepitaxial layer's free energy, three distinct 
growth modes can be observed. Technologically is most useful if the film grows in a layer-by- 
layer mode (Frank- Van der Merwe), when well controlled planar morphology is obtained. 
The deposited atoms diffuse on the surface and stick to the edge of nucleated islands. As 
a result the islands grow, finally covering the whole surface and completing the layer. On 
the top of the completed layer new islands start to nucleate, and the previous process is 
repeated. The growth has an oscillatory character in time, which can be observed using 
refiection high-energy electron diffraction (RHEED) or other experimental techniques. 

If the overlayer does not wet the surface, islanding is observed (Volmer- Weber mode), 
marked by the dumping of the intensity in RHEED measurements and absence of oscillations. 

And finally, if the overlayer wets the surface, but the overlayer stress is unfavorable, the 
film might grow in a layer-by- layer mode, followed by islanding (Stranski-Krastanov mode). 

One of the most widely studied heteroepitaxial structure is obtained by growth of Si on 
Ge, or Ge on Si. The Ge lattice is 4% larger than the Si lattice, thus generating considerable 
strain infiuencing the heteroepitaxial growth. As a result Ge grows on a Si (100) lattice in 
Stranski-Krastanov mode, while Si on Ge follows the Volmer- Weber model. 

The lattice mismatch generates islanding after 3-4 monolayers of layer-by-layer growth 
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during Ge deposition on a Si(lOO). Recently Copel, Renter, Kaxiras and Tromp proposed 
the use of a surfactant monolayer of As to reduce the surface free energies and effectively 
alter the growth mode . 

The microscopic mechanism responsible for the unusual effect of the As on the growth 
process is partially understood. The As layer, with one extra valence electron per surface 
atom, fills the dangling bonds of the Si(lOO) surface, creating a stable termination. Further- 
more, As segregates to the surface during growth. Si or Ge atoms deposited on the surface 
covered by an As monolayer rapidly exchange sites with the As and incorporate into the 
subsurface. As a result the heteroeptaxial structure incorporates negligible quantities of As. 

Two main mechanisms were proposed to explain the effect of the surfactants on the 
growth process |T5|JT7[1 . The first is a dynamic one, based on enhanced incorporation of the 
growth atoms. The As atoms drive any incoming Si or Ge atom to subsurface sites due to 
their ability to easily segregate. For a surfactant free surface the deposited atoms diffuse on 
the surface, until they reach a step or a defect, where they stick. In contrast, with surfactant 
the freshly arrived atoms are driven into subsurface sites by the exchange mechanism with 
the As atoms, their diffusion being severly curtained. Thus in the presence of a surfactant 
the adatom can be incorporated without a step or a defect. 

The second mechanism is using first principle calculations to explain the effect of the As 
atoms on the stress distribution of the surface layers. A shortcoming of this mechanism is 
that it predicts only an increase in the epitaxial thickness before islanding appears, but does 
not account for the change in the growth mode. Experimental results indicate that after 
about 50 ML the strain is fully relieved, supressing the driving force for island formation. 

Although probably the combination of the two effects are responsible for the supression of 
islanding, the nonequilibrium theory proposed in this paper is based on the first mechanism, 
providing a quantitative formulation of the dynamic phenomena occuring during surfactant 
mediated growth. 

In addition to the mentioned investigations [|T^|T3 5 ^ number of subsequent experiments 
showed that surfactants can change the surface morphology in a wide variety of systems. It 



was found that both Sb and As can act efficiently as a surfactant for the Si/Ge system [IT? 



The effect of the surfactant on the lattice strain and on the appearance of dislocations were 



studied in details p8|-pT[]. Low-energy electron microscopic observations were used to gain 
futher understanding in the local exchange mechanism between Ge and surfactant. It has 
been argued that surface energy anisotropy, instead of surface energy, is determining the 
changes in the growth mode of the Ge/Si system 

Further experimental investigations found that Sb alters the growth of Ag on Ag(lll) 
23| . Since submonolayer surfactant coverages were used, a new mechanism to explain these 



experiments was proposed. According to this the Sb attaches to the edge of the islands and 
lowers the interlayer diffusion barrier of the diffusing adatoms. The Sb is moving together 
with the edge of the growing islands and probably is segregated at the surface when the 
islands coalesce. 

In subsequent experiments antimony was found to change the structure of islands in 
Ge/Si growth and Te was used as surfactant to sustain layer-by-layer growth of InAs 
on GaAs(OOl) 

Since the most investigated system is the Ge/Si growth with As or Sb as surfactant, 
in what follows referring to the surfactant mediated growth we have in mind this system. 
Whether the proposed theory applies to all mentioned experiments, or additional effects has 
to be considered, is an open question, which will not be addressed here. 



III. NONEQUILIBRIUM THEORY 

In order to construct a nonequilibrium theory to account for the nontrivial effect of the 
surfactant on the growth, we have to study separately the dynamics of the interface and 
surfactant, and then consider the possible coupling between the two quantities. 

As mentioned above, under ideal MBE conditions, relaxation proceeds via surface diffu- 
sion. Atoms deposited on the surface diffuse until find an energetically favourable position 
(usually near a step or a dislocation), where they stick mostly irreversibly. The diffusive 
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dynamics is conservative, i.e. it does not change the mass of the film. The only change in 
the mass might come from deposition or desorption. 

In contrast to the ideal MBE, there is experimental evidence that surfactant mediated 
growth of Ge on Si proceeds by highly local Ge incorporation with minimum surface diffusion 
Ge atoms that adhere to the As-capped surface rapidly exchange sites with the As 



atoms and incorporate into subsurface sites. 

In the absence of surface diffusion, the growth equation may contain terms which violate 
mass conservation |^ . The simplest nonlinear growth equation with nonconserved dynamics 



was introduced by Kardar, Parisi, and Zhang (KPZ) |32| : 

dth = uV^h + X{Vhf + 7]. (1) 

Here h{x, t) is the height of the interface in d = d' + 1 dimensions. The first term on the 
right hand side describes relaxation of the surface by a surface tension u. The second term 
is the lowest order nonlinear term that can appear in the interface growth equation, and is 
related to lateral growth. T]{x,t) is a stochastic noise driving the growth; it can describe 
thermal and beam intensity fluctuations. 

Eq. ([I|) is the lowest order nonlinear equation compatible with the basic symmetries 
of a growing interface: it is isotropic in the substrate directions (x —x transformation 
leaves the system invariant), and invariant to translation both in the substrate directions 
(x x+a) and in the growth direction {h — > h+b). But there is a broken up-down symmetry 
in h: the transformation h ^ —h does not leave the system invariant. The explanation to 
this broken symmetry is based on the existence of a preferred growth direction for the 
interface. In the absence of the nonlinear term A this symmetry is obeyed as well. Another 
important property of this equation is that higher order nonlinear terms are irrelevant, i.e 
they do not effect the growth exponents (to be defined later). 

Additional terms in (0) will include the coupling to the surfactant fluctuations. 

In describing the dynamics of the surfactant we shall choose as parameter the width 
of the surfactant layer, v{x,t) (see Fig. 1). Throughout this paper is assumed that the 
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surfactant layer is very thin, thus nonlocal effects do not contribute to the dynamics. The 
typical experimental coverage, which is the spatial average of v{x,t), is around 1 ML. For 
coverages smaller than 1 ML holes might appear in the surfactant layer. Since the to be 
proposed growth equations do not depend in an explicit form of the thickness of the layer, 
but only on its spatial derivative, the system remains well defined even in the presence of 
such a holes. 

An efficient surfactant must fulfill two criteria: it must be sufficiently mobile to avoid 
incorporation, and it must surface segregate. Careful experimental studies showed for the 
Ge/Si system that the bulk As concentration is less than 1%; thus the effect of As on growth 
is a surface phenomena [lT5| , p!7| . 

Neglecting the desorption of the surfactant atoms, the equation governing the surfactant 
kinetics obeys mass conservation. 

This leads to the continuity equation 

dtV = -V-] + v\ (2) 

where t]' is a conserved uncorrelated noise which incorporates the random local fluctuations 
of the surfactant, and j is the particle-number current density. The simplest linear equation 
with conserved dynamics correctly incorporating the effect of surface diffusion is 

dtv = -KV\ + r/'. (3) 

Eq. d^) can be obtained from by using a current j ~ V/i, where /i is the local chemical 
potential on the interface. Considering /i ~ V^f , i.e. depends only on the local curvature 
of the thickness (describing local surfactant agglomerations), we obtain (j^). 

To account for the coupling between the growing surface and the surfactant it is necessary 
to introduce additional terms in Eq. (|I|) and (^. There are two main criterias which restrict 
our choice: The coupling terms must satisfy the symmetry conditions characteristic of the 
interface and the obtained set of equations should be self-consistent, i.e. the resulting 
dynamics should not generate further nonlinear terms. In addition the coupling terms 
included in Eq. (§) must obey the required mass conservation for the surfactant. 
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The simplest set of equations that satisfy the above conditions is 



(4a) 



dtv = -KV^v + iV^liVh) ■ (Vv)] + r/i, (4b) 

where the noise terms t]q and rji are assumed to be Gaussian distributed with zero mean 
and the following correlator: 

< r],{x, t)r],{x', t') >= V,5{x - x')S{t - t'). (5) 

Here 

Vo = Do (6) 

and 

Di = -DiV' + D2V\ (7) 
The D2 term is generated by Dq and Di as will be shown below. 

The generic nonlinear term (Vf)^ in ( Pa| ) can be derived using symmetry principles. In 
( ^bl) the V^[(V/i) ■ (Vt;)] term results from a current j = — V[(V/i) ■ (Vt;)], and obeys mass 
conservation. Geometrical interpretation of this term suggests that a positive 7 drives 
the surfactant to cover uniformly the irregularities of the surface, i.e. enhances the wetting 
properties [Q. A negative 7 has the opposite effect, assigning a non- wetting character to 
the surfactant. Since in experiments there is no evidence of surfactant agglomeration (non- 
wetting character), but it is energetically favorable to terminate the Ge layer with As atoms, 
we assume that the surfactant wets the surface, thus 7 > 0. 

The quantity of main interest is the dynamic scaling of the fluctuations characterized by 
the width [0] 

wl{t,L)= <[h{x,t)-h(t)f> =L'^-f(t/L^-) (8) 
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where xo is the roughness exponent for the interface h{x,t), and the dynamic exponent zo 
describes the scahng of the relaxation times with the system size L; h{t) is the mean height 
of the interface at time t and the <> denotes ensemble average. The scaling function / has 
the properties 

f{u ^ 0) ~ (9) 

and 

f{u—^ oo) ~ const. (10) 

In a similar way one can define Xi^nd Zi to characterize the fluctuations in the surfactant 
coverage f (x, t). 



IV. ANALYTICAL STUDY 

For /5 = 0, Eq. Pa| ) reduces to the KPZ equation For a one-dimensional interface 
the exponents can be obtained using DRG, resulting in the roughness exponent x = 1/2 and 
in the dynamic exponent 2; = 3/2. For higher dimensions unfortunatelly no exact results are 
available. But due to the non-renormalization of the nonlinear term A, the scaling relation 
X + z = 2 exists between the exponents, valid in any dimension. This reduces the number of 
independent exponents to one. A number of conjenctures exist in the literature regarding the 
higher dimensional exponents, but so far none of them is proved. But numerical simulations 
on discrete models and direct integration of (|I|) helped to obtain reliable estimates for the 
exponents in higher dimensions as well. For the physically relevant dimension, c? = 2 + 1, 



extensive numerical simulations give xo = 0.385±0.005 and zq ~ 1.6 |3^. Thus the interface 
is rough and the roughness increases with time as Wo(t) ~ txo/^o_ 

For 7 = 0, Eq. (|¥5D is the fourth order linear diffusion equation with conserved noise 
(|), which can be solved exactly, resulting in zi = 4 and Xi = p8| , |39| . In d = 2 + 1 these 
exponents do not change even if additonal nonlinear terms, compatible with the symmetries 
and conservation laws of (H), are added to the linear equation. 



Thus, neglecting the couphng terms, Eq. ( ^al) and ( pq) predict rather different values 
for Zi and the roughness exponents Xi- To see how the couplings change this behavior we 
have investigated Eq. using a DRG scheme. 

For this we rewrite Eq. (11) in its Fourier components 



h{k, oj) = fjQ{k, uj)Go{k, u) — XGo{k, ^) J J df'QdVt q{k — q)h{q, Q)h{k — q,uj — Q) 

-(3G{k,uj) f f d'^qdn q{k-q)v{q,n)d{k-q,uj-n) (11) 



d{k, u) = fii{k, uj)Gi{k, u) + -fk'^Gi{k, u) J J d'^qdQ q{k - q)h{q, n)v{k -q,u-n) (12) 

where fji{k,uj), h{k,uj), and v{k,uj) are the Fourier components of the corresponding 
quantities and the correlators have the form: 

During the DRG calculations only one dynamic exponent z = zq = Zi was used, valid if 
the equations (^) do not decouple. Equations (|TI],|T2p are the starting point for the pertur- 
bative evaluation of h{k,uj) and v{k,uj). The basic diagrams are indicated in Fig. 2. The 
fast modes are integrated out in the momentum shell e~^Ao < \k\ < Aq, and the variables 
are rescaled as x ^ e'x, t e^H, h e^°''h, and v — e^^'f . The calculations have been 
performed up to one-loop order. 

In what follows we shall skip most of the details of the calculation, the interested reader 
is referred to the literature We shall presents only the main parts which are relevant 



to further arguments. 

The first result is that the diagrams contributing to A cancel each other, resulting in the 
flow equation 

^ = X[z + xo-2] (15) 
10 



providing us with the scahng relation 

2; + Xo = 2. (16) 

This relation is known to be the property of the KPZ equation and it is a consequence of 
Galilean invariance (GI). Since the DRG conserves the GI, this scaling law is expected to 
remain valid to all orders of the perturbation theory. 

A second scaling relation can be obtained from the non-renormalization of the diffusion 
coefficient D\ : 

'^ = D,[z-d' ~2~2x,i (17) 

resulting in 

z - 2xi - rf' - 2 = 0. (18) 

The diagrams that contribute to Di (see Fig. 3) have a prefactor proportional to k^, thus 
they are irrelevant (k is the wave vector in the Fourier space). They in fact contribute to 
D2, justifying its introduction in (|^. 

These two scaling relations already indicate that the coupled interface/ surfactant system 
is qualitatively different from the uncoupled one. For a planar interface {d' = 2) (|l^,|l^) give 

Xo + 2xi = -2, (19) 

which means that at least one of the exponents has to be negative. 

A third scaling relation unfortunately is not available, but insight can be obtained from 
numerical integration of the flow equations obtained from the DRG. A correct flow must 



not scale the nonlinear terms (3 and 7 to zero, which would decouple Eq. ( ^a]) and (^). 
The finiteness of the nonlinear terms guarantee the validity of the scaling relations (|I6|,|T8|) 
as well. The integration showed the existence of two main regimes: 

(i) In the first regime one or both of the coupling terms {(3, 7) scale to zero. In this case 
the two equations become completely (both coupling terms vanish) or partially (only one 
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coupling term vanishes) decoupled, and the two equations might support different dynamic 
exponents z. The DRG scheme used is not reliable in this regime. 

(ii) The presence of a strong coupling fixed point is expected when both of the nonlinear 
terms diverge. The integration shows that this coupled phase exists only for z > 8/3. The 
coupled phase is stable against small fluctuations in the coefficients and exists in a finite 
region of the parameter space. Since under experimental conditions small fiuctuations in 
the value of the control parameters are always expected, the stability of the system against 
them ensures the persistence of the coupled phase. But large deviations of the parameters 
introduce instabilities, which result in the breakdown of the smooth phase. This is in accord 
with the experimental observation, that surfactant induced layer-by- layer growth develops 
only under well controlled experimental conditions. 

It is important to note that although there is no identifiable fixed point, in this phase 
the scaling relations ([1^, |1^) are exact. According to ([T6| ) for z > 8/3 the roughness 



exponent of the interface Xo is negative (see Fig. 4). With a negative roughness exponent, 
every noise-created irregularity is smoothed out by the growth dynamics and the resulting 
surface becomes fiat. Thus the coupling of the surfactant to the growing interface results 
in the suppression of the surface roughness. This corresponds exactly to the experimentally 
observed behavior, i.e. the addition of the surfactant suppresses islanding, resulting in a 
morphological transition from rough (without surfactant) to fiat (with surfactant) interface. 

The roughness exponent of the surfactant from Xi (HH) is negative if 2; < 4, while for 
2; > 4 it becomes positive (See Fig. 4). In the Ge/Si system, for example, the As has a 
saturation coverage of 1 ML, which is independent of the system size and is governed only 
by the microscopic bonding of the As to the Ge dangling bonds. One expects no relevant 
fiuctuations in the thickness of the coverage; this requires a negative roughness exponent for 
the surfactant and thus limits the dynamic exponent to values smaller than four. 

The DRG analysis fails to provide the exact value of the dynamic exponent z. As in 
the case of many other growth phenomena, simple discrete models might be very helpful to 
obtain its value (see discussion later). Summarizing the results of the direct integration of 
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the DRG equations, for z > 8/3 the existence of a strong couphng fixed point is observed, 
in which the interface roughness exponent is negative, corresponding to a flat phase. There 
is no upper bound in z for the existence of this phase, but physical considerations suggest 
that z < 4, in order to allow the uniform surfactant coverage observed experimentally. 

V. CONCLUSIONS AND FUTHER DEVELOPMENTS 

In the previous sections we introduced a set of coupled equations compatible with the 
basic symmetries and conservation laws of the surfactant /interface system studied exper- 
imentally. The main feature of these equations is that they predict a negative roughness 
exponent. We have argued that a negative roughness exponent describes a flat interface, in 
accord with the experimental observations. A natural question arises here: Is there any pre- 
dicting power in this theory, or just reproduces the experimental results without generating 
further inquiries? 

In this section we examine the predictions made by the theory. The limits are presented 
as well: what are the physical ingradients we neglected, and whether and how could they 
be incorporated in a new theory along the presented lines. 

As we have noted earlier, the analytic study does not provide us with the exact value of 
the exponents. But predicts that the dynamic exponent z lies in the narrow range between 
8/3 and 4. If we could measure somehow the dynamic exponent z, the scaling relations would 
provide us with the other exponents. In fact, if one would be able to measure experimentally 
any of the exponents Xi? A? or z, the other exponents could be obtained via ( p!6| , |l8D . 

The scaling theory (|,|],|10D predicts that an originally flat interface becomes rough as a 
power law of time, w ~ t^. Since in our case (3 is negative, an originally rough interface 
becomes smooth as a power law of time, until a limiting small roughness is reached. The 
only difference is in the system size dependence of the roughness: while in the usual growth 
models the roughness increases as a power of L, in our case the interface is smooth, with a 
small thermal roughness Wq, independent of the system size. Thus a possible experimental 
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check of the previous predictions would start from an initially rough interface and monitor 
directly the decrease of the roughness in time and fit the obtained curve with a power law. 
Previous experimental results indicated that it is possible to obtain the time dependence of 
quantities directly related to the surface roughness []5T|-[55|]. It would be interesting to see 
whether for the surfactant system such a study could be carried out. 

Such an experiment would result in the exponent (3o for the interface (/?o = Xo/^): from 
which using the scaling relations z and xo could be determined. Hopefully the determined 
z would fall between the boundaries predicted by the theory. 

Futher test of the theory might come form the direct numerical integration of the coupled 
equations (Q), with the aim to look for the coupled phase and obtain the value of the critical 
exponents. Integration proved to be successfull in obtaining the exponents for the KPZ 
equation , and for checking the DRG results for other coupled systems . 



Constructing and investigating discrete models in the same universality class as the 
studied continuum equations is another efficient and frequently very accurate way to obtain 
the scaling exponents PS|-BD],P3| . For nonconserved coupled equations (see later) such models 



have been investigated and gave results in accord with the DRG and numerical 
integration |^ . 



And finally let us mention some open problems related to the presented theory. It is 
important to note that introducing Eq. (H) we did not use directly the existence of the 



strain which appears due to the lattice mismatch. Although an important problem a 
continuum description of strain-induced roughening is still missing. The proposed model 
is expected to describe the con^^/ec? surfactant/interface system, but decoupling the surfactant 
does not necessary result in an equation describing heteroepitaxial islanding. Further studies 
are necessary to understand the microscopic (perhaps strain induced) origin of the nonlinear 
coupling terms. 

In Eq. d^) the desorption of the surfactant atoms is neglected by considering that (|4ED 
obeys mass conservation. Lifting the conservation law, (^^ should be replaced by a non- 
conservative equation. Such a system has been recently studied [^ , |47|| , and it was found 
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that in most cases the couphng does not change the KPZ scahng exponents. Enhancement 
of the exponents is possible only when the couphng is one-way, i.e. one of the equations 
decoupled from the other one is acting as source of correlated noise. 

Further linear and/or nonlinear terms added to (^) might influence the dynamics of the 
system. The goal here was to derive the simplest set of equations predicting the experimen- 
tally observed morphological phase transition; the study of other possible nonlinear terms 
and their relevance is left for future work. 

Another shortcoming of the presented theory is that it does not predict oscillations in the 
interface roughness in the layer-by-layer growth regime, as is expected experimentally. This 
is due to the fact that the present continuum theory does not account for the discretness of 
the lattice, responsible for the oscillations. But such a discrete pinning potential in principle 
can be introduced in (^. The effect of such a pinning potential for both the conserved and 
nonconserved equation was studied in the literature [p7| , |58| . It would be interesting to see 
how the coupling terms interact with the lattice potential, and whether such a calculation 
leads to a coupled phase with periodic oscillations in time. 

In conclusion, I have introduced a set of equations to describe the interaction of a growing 
surface with a surfactant. The main experimentally motivated requirements for were: (a) 
no surface diffusion of the newly landed adatoms; (b) conservative and diffusive surfactant 
dynamics, originating from neglecting incorporation and desorption of the surfactant during 
the growth process. The obtained equations indicate the existence of a coupled phase, in 
which two scaling relations between the three exponents are available. In this phase, the 
roughness exponent of the interface is negative, morphologically corresponding to a fiat 
interface, as observed experimentally. 

Moreover, Eq. (|^) serve as a good starting point for future studies of an interface coupled 
to a local conservative field, a problem of major interest in the context of recent efforts to 
understand the general properties of nonequilibrium stochastic systems. 

I wish to thank M. Gyure and E. Kaxiras for useful discussions and comments on the 
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FIGURES 

FIG. 1. Schematic illustration of the studied surfactant/surface system. The figure represents 
a cross section of the two dimensional surface of heigh h{x,t) covered by a thin surfactant layer 
with thickness v{x,t). A newly arriving atom penetrates the surfactant and is deposited on the 
top of the growing interface h{x,t). 

FIG. 2. Diagrammatic representation of the nonlinear integral equations (11,12). 

FIG. 3. The leading contribution to the effective noise spectral function. The encircled noise 
term corresponds to X>i, and consists of Di and D2 according to (7). 

FIG. 4. The dependence of the roughness exponents xo and xi on the dynamic exponent z, 
according to the scaling relations (16) and (18). 
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